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1 Introduction 



D-branes in bosonic string theory are unstable, and their decay provides a relatively 
clean non-trivial example of a time-dependent background in string theory. In super- 
symmetric string theory, there are also unstable D-branes. The stable D-branes carry 
charges, so a pair of oppositely charged branes at subcritical separation also becomes 
unstable and will decay in a similar manner. Investigations of properties of stable 
D-branes have lead to many important insights, so one may hope that the same will 
be true for unstable D-branes. The annihilation of an unstable pair of D-branes is a 
basic process to be understood, and also a building block for many string-based cos- 
mological scenarios. However, D-brane decay is computationally a more demanding 
process to study. Many simple questions remain unsolved, and have received little 
attention recently^ In this paper we investigate the amplitude for the emission of 
a pair of closed strings from a decaying brane, in the simple "half S-brane" rolling 
tachyon decay background [4,5]. We focus on Dp-branes with p < 25. 

The single closed string emission channel is relatively well understood, both in 
bosonic [6-8] and superstring [9] theory, also in the presence of background electric 
fields, fluxes, or winding [10-13]. Bulk-boundary amplitudes (emission of a closed 
string from a perturbed D-brane) can also be calculated analytically in bosonic and 
super symmetric theories [14-16] (see also [17,18]). 

For multi-string amplitudes, one strategy to calculate the associated n-point (bulk 
or boundary) amplitudes was based on previous success in bulk Liouville theory [19]. 
Also in boundary Liouville theory, two-point and three-point functions have well- 
defined analytical expressions [20-22]. In [23,24] it was proposed that analytic con- 
tinuation in coupling b — > i and field space <fi — > iX° to timelike Liouville theory 
could be used to obtain results in the rolling tachyon background. For bulk-boundary 
amplitudes, analytic continuation from spacelike theory lead to a result [15,17] which 
is in agreement to that of [14] which was based on a more straightforward approach. 

Closed string pair production, or associated two-point functions were investigated 
in [23, 24] but the obtained amplitudes did not display the expected pole structure. 
Ref. [14] analyzed the problem and identified some potential problems in the ana- 
lytic continuation. For additional related work, see [25]. We will comment on the 
analytic properties in more detail in Section [2j For the "full S-brane" background, 
an interesting prescription to compute disk n-point amplitudes was proposed in [26], 
and its higher order extensions were investigated in [27]. It would be interesting to 
investigate this method further to compare the results with those obtained on the 
half S-brane background. Additional work on string scattering from decaying branes 
can also be found in [28]. See also the reviews on open string tachyon dynamics [29] 

1 Except for the more ambitious open string field theory program, where notable progress in 
tachyon condensation was obtained in [1,2] and subsequent work, see [3] for recent discussion in 
relating the results to boundary conformal theory. 
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Figure 1: Two channels for emission of a pair of closed strings from a decaying brane. 
Note the t <-> s difference compared to string scattering, see Fig. 2 of [37]. 

and on Liouville field theory [30]. 

In this paper we continue developing the approach initiated in [14], in particular 
we use the relation of correlation functions in the rolling tachyon background to 
thermal expectation values in a classical log gas of unit charges in two dimensions. 
The thermal interpretation was further explored in [31-36]. We also use the powerful 
contour integral representation which was developed in [25] to study n-point boundary 
functions. Our main result is an amplitude (obtained by working in the high-energy 
limit) for the emission of a pair of closed string tachyons from the decaying brane. 
The formula differs from the previously obtained ones by its more promising pole 
structure. It contains the two expected channels: an s-channel, where an on-shell 
closed string propagates first before splitting to the two emitted closed strings, and 
a t-channel, where a pair of on-shell open strings propagate first before closing off to 
the pair of closed strings, depicted in Figure [TJ 

Our new calculational ingredient is the use of electrostatics to estimate the leading 
contribution to the correlation function. While electrostatic analogues have a long 
history in string theory [38] , our approach is reminiscent of that by Gross and Mende 
[39] (see also [40,41]), who used electrostatic equilibrium conditions to identify the 
dominant saddle point contributions to string scattering amplitudes at arbitrary order 
to determine their high energy behavior. We will also be interested in the high energy 
limit (for the emitted strings), but we consider only the leading disk amplitude, which 
is already complicated due to the rolling tachyon background which will correspond to 
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a thermal charged gas background. In the high energy limit, our amplitude displays 
the expected exponential ~ e~ 2nuJ suppression, with the sum of the energies of the 
emitted strings in the exponent. 

The paper is organized as follows. Section [2] contains a brief derivation of the 
perturbative series expansion of the bulk two-point amplitude. We then comment 
briefly on previous work and analytic properties of the series expansion, and convert 
it to an integral representation appealing to [25]. Section [3] presents the electrostatics 
approach and applies it to the emission amplitude calculation. Section [4] introduces 
Mandelstam variables and presents the result for the amplitude formula and identifies 
the poles. We then consider its high energy asymptotics in different kinematic limits 
and find that the amplitude becomes at least exponentially suppressed, as expected. 
In the end there are three appendices, Appendix [A] collects some details on kinematics 
and the amplitude calculation, for readers' convenience. Appendix [B] discusses an 
alternative attempt to obtain an exact result for the two-point amplitude, using a 
Schur polynomial method. Appendix [O contains some details of the electrostatics 
calculations. We also plan to present a more general discussion of the electrostatic 
approximation in a forthcoming work [42]. 



2 The bulk two-point amplitude 

We begin with setting up the notations for different correlation functions. We work 
in units where a' — 1. The correlation functions are worldsheet path integrals 

A 2 (K,ULi) = / fl^(flV(w a ,w a )e- ss ^\ , (1) 

a=l \a=l / f rcc 

where the boundary deformation representing the half S-brane [4, 5] is 

6S hdly = X j dte x ° {t \ (2) 

The vertex operators for the closed string tachyons are 

V{w a ,w a ) =e lka - x{Wa '™ a) , (3) 

where a = 1,2. We adopt a notation £ a = — iu a and break up the spatial momentum 
to parallel and perpendicular directions to the unstable Dp-brane: k a = (hi, k^). We 
also denote kl = (u a , ki) = (z£ a , kl). On-shell conditions for the bosonic closed string 
tachyons are kl = ^ + (k a ) 2 = 4. Furthermore, one finds that the overall (spatial) 
parallel momentum is conserved: kj + k 2 = 0. 

The worldsheet correlation functions can be evaluated by first isolating the zero 
modes from the oscillators, X M = + X' M , and then expanding the boundary defor- 



4 



mation into a power series. This yields 



N ' 2 d 2 w n 



A 2 (K,fcJ) = dx^e^^^^m /n 

AT=0 ' 0=1 

N , / N 2 

x JJ — ( JJ e x '° {tl) JJ e**--^' 

i=l \ i=l a=l 



2tt 



(4) 



where we introduced z = 27c\e x . 

After working out the contractions, with the Neumann and Dirichlet Green's 
functions (see Appendix lAl and, e.g., [14] for more details), (BJ becomes 



A 2 ({uo a ,k a }) 
A 2 (x°) 



dx°d p x ll e i ^ k « x »A 2 (x° 
r 2 j2„.. 

/ n 



(5) 



N=0 



0=1 



2tt 



^11- WaiSil J^-^*>Z a ({u; a) U; N) 



(6) 



a,b=l 



where 



AH / AA 2?r 

1=1 l<i<j<N 



- w a e 



(7) 



0=1 



2.1 Some general comments 

We fix the residual global conformal symmetry on the disk by placing the other 
bulk vertex operator to the origin (w 2 = 0) and the other one to the real axis at 
1 > W\ = r > 0. Let the (imaginary) energy of the closed string vertex at W\ — r be 
£i = £ and the energy of the string at w 2 = be £ 2 - The first challenge in calculating 
the amplitude is the integrals (j7|), which now become 



Z^;N) = - 



N 



.i=l 



dU 
2?r 1 



n 

X<i<j<N 



We can make the following general remarks: 



For < r < 1 the integral Z 2 (r, £; N) is an analytic function of £. 

For negative integer £ and for iV > |£| the integral can be evaluated by using 
identities of the Schur polynomials [24] giving Z 2 {r, i;N) = (1 - r 2 )-« . For 
positive integer £ the integral is a polynomial in r, and may be written down 
explicitly by using Schur polynomials (see Appendix IB]) . 
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• In the limit N — > oo authors of [14] also find Z 2 (r, £; AT) — > (1 — r 2 ) ^ 

• For r — ► 1, and if Re£ > —1/20 Z 2 (r, £; iV) approaches smoothly the boundary- 
one-point and bulk-boundary functions [14] (see also [15]). 

In the final amplitude the interesting region is £ = — za;, where u is the energy of 
the closed string at w\ = r. In the case of (integrated) boundary amplitudes, it is 
possible to analytically continue from positive Re £ to imaginary values [25] . The bulk 
integrals approach smoothly the boundary ones for positive Re£ in the limit r — > 1, 
so it is natural to assume that the Re£ > region may be a good starting point to 
study the complete bulk amplitudes. 

There are closely related problems where the scattering amplitude, the sum of iV 
over an infinite series of terms, is essentially equal to a generic iVth term after analyt- 
ically continuing iV — > —i ^2i a ^ a - in spacelike Liouville theory, for bulk amplitudes a 
related observation was made in [43]. In [14] it was observed for the bulk-boundary 
amplitude in the decaying D-brane background. Then, ref. [25] constructed a detailed 
contour integration method to show how boundary amplitudes reduce to the analyt- 
ical continuation of a term. Following [14, 25] we expect that the bulk amplitude is 
given by a certain analytic continuation to negative N: 

A 2 = "^T^ / 1 rfrrr^^l-r 2 )K^)M^) Z2(r ^. iV = ^_ 6) . (g) 
sm7r(£ + £ 2 J Jo 

In the integrand of (JHJ), at this point we only know the term Z 2 through its 
definition in (jSJ), which is quite complicated. However, it allows us to make some 
observations about the pole structure of the amplitude. There singularities at each 
endpoint r = 0, 1 in the r integral give rise to pole structures, as in scattering from a 
stable D-brane [44], see [37] for a review J§ From (jSj) we can see that Z 2 is real analytic 
at r = and even in r. It is also analytic at r = 1 for Re£ > —1/2. Consequently, 

the singularities at the endpoints originate from the factor rr kl ' k2 (1 — r 2 ) 2 ^ ^ ^ ) , 
which is the same as in scattering from the stable D-brane. Let us consider the poles 
from r = in more detail. We may expand Z 2 as a Taylor series at r = 0, 

Z 2 (r,£;N=-£-6) = ^ 2 (0,£;iV=-£-£ 2 ) + ir 2 |^Z 2 (r,£;iV = -£-£ 2 )| r=0 +--- 

= l + a 2 (£,6)r 2 + --- , (10) 



2 Authors of [14] use functional relations which hold only for Re£ > —1/2 in their discussion of 
the r — ► 1 limit. 

3 On the other hand, scattering from a stable D-brane corresponds to the N = term in the series 
expansion (TJ]). 
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where a 2 (^,^ 2 ) is the proper analytic continuation of the series coefficient^ to iV 
— £ — £ 2 . Integrating over r term by term in (jH]) then yields 



7r(27rA)- 



sin7r(£ + £ 2 ) 
sin7r(£ + £ 2 ) 



1 

drrr kl ' k2 
1 

h-k 2 + 2 



+ a 2 (£,6)-|(M' 



fci • k 2 + 4 



'11} 



The poles are located at k\ ■ k 2 = —2, —4, . . ., as in scattering from ordinary D- 
branes [37, 44] . The above calculation highlights the structure of the integrand at 
r = 0. 

Similarly one can check that the singularity at r = 1 gives rise to poles at 2{k\) 2 = 
2, 0, —2, . . . where we used the mass-shell condition k\ = 4. 

The remaining paper focuses on developing an electrostatic approximation for Z 2 . 
It is analytic at r = 0,1, so the poles indeed arise as described above, and we can 
verify the expressions explicitly. We shall discuss in detail only the case £ = £ 2 where 
the production amplitude is symmetric^ 



3 Electrostatic calculation 

Our starting point is to view [14,16,25,31] Z 2 as the partition function of a Coulomb 
gas [45], with two additional external charges £, £ 2 at fixed positions. We write it in 
the form 

mUhN) = ±jfl d ±e-™, (12) 
where the Hamiltonian H reads 

N 2 

H=- £ log |e* - e*i I - ]T £a |e* - , (13) 

l<i<j<N i=l a=l 

and we fixed the inverse temperature to /3 — 2. As seen from the Hamiltonian, we 
have an ensemble of iV unit charges on the circle at e ltl , . . . , e ltN , which interact via 

4 See Appendix [5] for a calculation of the series coefficients for N = 0, 1, 2, . . .. In particular, for 
integer N one can extract from eqn. (|7ip that 02 = £ 2( 5jv>o- 

5 The reason for this restriction is a puzzling observation for which we do not have a good ex- 
planation at the moment. When we use conformal symmetry to place the other vertex operator to 
origin and the other one to w = r, the amplitude ©-([7]) becomes asymmetric with respect to the 
exchange 1 <-> 2 of the strings, if they have different energies, £ ^ £2- As a result, the poles seem not 
to respect the exchange symmetry. So, the use of conformal symmetry to fix the vertex positions 
in the rolling tachyon background involves a subtlety, unless £ = £2, in which case the strings are 
exchange symmetric "in the time direction" . 
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the Coulomb potential^ They are subject to an external electric field created by two 
fixed charges. 

A useful strategy is to study the system using classical electrostatics, first in 
the limit of large number iV of unit charges. In this limit the unit charges can be 
treated as a continuous charge distribution. Its density is modified by the external 
charges £ a and we are interested in finding the electrostatic equilibrium configuration. 
Notice that the present system is analogous to the circular unitary ensemble [5] of 
random matrix theory, where the connection to the Coulomb gas is a standard tool. 
Our method indeed resembles the electrostatic derivation of Wigner's semi-circle law 
[46] (see also [47]). It is also reminiscent of the work of Gross and Mende [39], 
who used electrostatic equilibrium conditions to identify the dominant saddle point 
contributions to string scattering amplitudes at arbitrary order to determine their 
high energy behavior. We will also consider the high energy limit (for the emitted 
strings), but we restrict to the leading disk amplitude, which is already complicated 
due to the charge density background. A more detailed description of our method 
and other applications will be discussed elsewhere [42], here we focus on aspects that 
are relevant for the bulk two-point amplitude. 

We use a saddle point approach and write 

log Z 2 ({Q; N) ~ -2H ~ -2S tot , (14) 

where H Q is the minimum value of the Hamiltonian (flBl and £ tot is the (total) electro- 
static energy of the corresponding equilibrium configuration with continuous charge 
densities. The string production amplitude (jHJ) requires N = — £ — £ 2 , so that we 
need to study the electrostatic problem in the regime N ~ £ ~ £ 2 - It turns out that 
nontrivial electrostatic configurations are indeed obtained when the external charges 
are of the same order as the sum of the unit charges. Slightly more precisely, the 
electrostatic configurations of interest to us are in the limit iV — > oo with the ratios 
N/£ and N/& fixed0 

Let us then discuss how the total electrostatic energy £ tot is calculated. The first 
ingredient is the (complex) potential V(w), the physical potential is the real part 

U(w) = ReV(w). (15) 

The potential is sourced by the continuous charge density p(4>) (by which we approxi- 
mate the N unit charges in the iV — > oo limit, with the normalization N = § d<pp((f))) 
and the external charges £ a , 

V(w) = - J d<pp(<j)) log (w - e^) - £a log(w - w a ) + const. . (16) 



6 Other unstable D-brane backgrounds correspond to different Coulomb gases [32-34]. 
7 Recall that the behavior of Z2 in the limit N — > 00 with £ a and w a fixed can be calculated by 
using the asymptotics of Toeplitz determinants [14, 48-50] . 
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We choose the zero level of the potential by setting the last constant term to zero. 
We then need to find V(w) and p((p) in the equilibrium configuration. We are inter- 
ested in the configurations where pU)) is strictly positive within a domain C of the 
unit circle and vanishes elsewherejfl Then solving the equilibrium configuration re- 
verts to a standard potential problem of classical electrostatics including the external 
charges £ a and a conductor which fills the domain C. In particular, in the equilibrium 
configuration the physical potential U(e 1 ^) must take a constant value, 

t/(e # ) = U , (17) 

in C: If the potential would vary within C, the charges on the conductor would be 
subject to forces and move until the potential adjusts to a constant value. 

The electrostatic energy is defined as the energy of the system in the equilibrium 
configuration, and may be expressed in terms of po(4>) as 

£tot = -^y#i# 2 po(0i)po(0 2 )log|e^-e^| 
/ d#o(0)log|e^- Wa | 

= £ c + £ c ^ . (18) 

We again set a possible constant contribution to zero. Then the energy matches with 
the continuum limit of the Hamiltonian ([13]) and thus approximates the saddle point 
value of — |logZ 2 . By using the fact that the physical potential £/(e 1 *) is constant 
within the domain C of nonzero po, we can express the result as 

N 1 

Stot = -Uq + -£c_ { , (19) 

where the first term includes one half of the "interaction energy" £ c _£. Hence, to 
calculate £ tot and the leading behavior of Z 2 for large N ~ £ ~ £ 2 , it is sufficient to 
find Uo and £ c _£. We present a general method for this in Appendix O 

In the next two subsections we work out the electrostatics approximation of the 
(two-point) bulk amplitude explicitly. We can omit the charge £2 at the origin since 
it does not affect the calculation. The treatment of the other (real) charge £ at the 
real line with < w\ = r < 1 is divided to two regimes, depending on the values 
of r and £/N. For small values of r the charge £ affects the charge distribution on 
the unit circle less, and it remains strictly positive for all <fi = argw. For r — ► 1 
(assuming £ > 0) the charge creates a gap in the charge distribution on the circle 
in the interval — d> r < d> < d) r n We discuss first the small r case. 



3 In this paper it is sufficient to assume that the domain C of non-zero p(4>) is connected, that is, 



an arc 
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The distribution remains strictly positive if £ is negative and N > |£|. We shall consider the 
gap creation only for positive £ where the r — ► f limit of Z2 matches with the one-point boundary 
function [14]. 
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3.1 Small r 



For small r one needs to solve the potential problem on the unit disk (without gaps), 
which is easily done by using a mirror charge. One can check that U (e^) = Uq = in 
this case when the constant term in (fT6l) is set to zero. At electrostatic equilibrium, 
the potential becomes 

V(w) = -£log-^^. (20) 
1 — wr 

The first part of (fl6|) . the contribution from the continuous distribution on the circle, 
becomes V c (w) = +£ log(l — wr) at equilibrium, which in turn equals the potential 
of the mirror charge at w = 1/r. As mentioned above, we dropped the trivial con- 
tribution — ^logw due to the charge at the origin which neither affects the charge 
densities nor the total energy. The equilibrium charge distribution on the unit circle 
is 

pm = i (n + a - jl r wu,. = i (s + { - 1 - ; :; 2 cos J , pd 

where the first (constant) term is fixed by § d(f)p(<f)) = N and the second (0 dependent) 
term is the charge density induced by the charge £ at W\ — r. The absence of the gap 
requires that p is positive, i.e., that 

N 

r < = r c , (22) 

which defines what we mean by "small r" . The total energy at equilibrium becomes 
£tot = 1 -S c ^ = ^V c (w = r) = ^log{l-r 2 ) . (23) 

Notice that in the large A^ limit with fixed £ we have r c — > 1, and the result ( 1231) is 
valid for all r. Hence in this limit we may use the result directly in ([141) , 

Z 2 (£- N) oc e~ 2£tot = (1 - r 2 Y e ; N -> oo , (24) 

which indeed reproduces exactly the large A^ result of eqn. (5.13) in [14]. However, 
as we will see next, at large r we obtain in general a different result. 



3.2 Large r 

As we saw in ( |22|) . for r > N/(N + 2£) there is a gap on the unit circle charge 
distribution. The electrostatic problem becomes more difficult. We present a general 
method to find the equilibrium charge density in the presence of gaps in Appendix [Cl 
Here we focus on the results in the case of one gap n = 1 and one external charge 
n — 1 which is relevant at large r. 
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Appendix O performs the calculation in a simplified geometry, where the unit disk 
(in w coordinates) is mapped onto the upper half plane (z coordinates) by 

1 — w 1 — w 

z = q (w) = = i cot (0 c /2) — — , (25) 

1 + w 1 + w 

where <p c will parameterize the location of the gap in the charge density. In the 
electrostatic problem (on the w plane), there is an induced image charge —N — £ at 
infinity, so that the net charge is zero. Following the point at infinity in the mapping 
q in (1251) . the image appears as an additional charge at z — q(oo) on the z plane. In 
the notation of Appendix O the actual charge is £ n = £i = £ and the image charge is 

£1 = £ at z — i8(r)/3 (w = r) 
£2 = —N — £ at z = -i(3 (w = 00) , (26) 

where (3 = cot(^ c /2) and 

8{r) = . (27) 

1 + r 

The asymptotic behavior of the electric field leads to a constraint (eqn. fl84"j) in 
Appendix [C]) which gives a relation 

or equivalently 

cos 2 ^ = l ~ X " = (29) 
2 l-tf(r)> i_(l_r) 2 ' (29) 

where \ = £/(N + £). Notice that <p c vanishes at r = r c = N/ (N + 2£) which signals 
the disappearance of the gap. By using eqs. and (19"01 we find for the constant 
value of the potential 

U = -—log- — — log—-—, 30 

2 1 — o[r) 2 1 + d(r) 



while the interaction energy of (l9Tj) is given by 

Be-* = aN+o 



log — — h x log ■ 



i + x & (i + x)(i + W) 



(31) 



The total energy thus reads 

£tot = + 2^ c ~ 5 



(iV + 20 2 , l + x ^ 2 l~x J\ 4 X 



4 to l + 5(r) 4 to l-5(r) 2 to (1 + 5{r)) 
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4 High energy behavior 



Let us then apply the electrostatic results to estimate the string production amplitude 
of (jHJ). In the bulk case there is a subtlety with analytic continuation of [14,25] 
N — > — £ — £ 2 - When r is kept fixed, the integral can be approximated by the "large 
r" result (1321) for N < N c , where N c = 2£r/(l — r), whereas "small r" approximation 
(I2"3"j) works for iV > N c . Thus there is a discontinuity at iV = iV c , which leads to 
apparent problems with the analytic continuation. The sum that we should calculate 
reads 



oo oo 
N=0 N=0 

,.o 



A 2 (z) ~ V(-^Z 2 (r,e;JV) ~ V(-z)^e" 2 ^ , (33) 



where z = 27r\e x and £ tot is given by either (J32D if iV < N c or (J22D if N > N c . 
The sum is well convergent and insensitive to the actual value of N c for \z\ < 1. For 
continuation to |z| > 1 we would like to apply the contour integral method of [25]. 
However, the discontinuity at N = N c seems to be problematic, if r is held fixed. 
The sum A 2 (z) is sensitive to N c and consequently diverges fast for z — > oo with an 
arbitrary proportionality constant, whence it is not possible to integrate the amplitude 
over x°. Working directly at imaginary £ = — iu might improve the convergence, but 
the analytic continuation is still tricky. 

Remarkably, the problems disappear if r is integrated over first. Hence we define 
the integrated partition function 

1 rr c ,1 / j , \ (N+20 2 /2 

drk(r)Z 2 (r,£;N) ~ / drk(r) (l - r 2 ) 4 + / rfrfc(r) ' v 



1 + «J(r) 



1-X \ / 4 X 



l-5(r)y V(i + ^(^)) 2 
= ^(C;iV), (34) 

where r c = N/(N + 2£) and fc(r) can be any arbitrary analytic function, but we will 
choose it appropriately below in (1361) . Now Z is an analytic function of N, and can 
be continued analytically to all complex iV plane. It has a few branch cuts: branch 
points are found whenever the endpoint of the integration, r c , hits singularities of the 
integrand. Branches are given by the various winding possibilities of the integration 
path. The singularities are found at r = 0, ±1, and oo, giving for the branch points 
N = 0, -2f , and oo. 

It is straightforward to check that Z grows only as a power law as iV — * oo in 
any direction on the complex plane. Hence, following the analysis of [25], we can 
apply the contour integration method, which gives our conjecture for the high energy 
behavior of the bulk two-point function 



A 2 c f dx° e^>° jT(-z) N Z(N) = m N=-t-&), (35) 
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where we fixed (see (Q) 



k(r) = rr kl ' k2 (l - r 



(36) 



and omitted the trivial factor 5 [k^ + fc 2 J • Explicitly. 



An 



-C-6 



drr 



sin7r(£ + £ 2 , 
: / rfrr fcl ' fc2+1 (1-r 2 ) 



(1 



r 2 ) 2 " 



^2 +(fc ll ) 2_ (fc X ): 



(37) 



2\M( fe i) 2 -( fe i x ) 2 



X 



(l + ^(r)) s 
4x 



e 2 



1 + <5(r) 



(C-6) 2 /2 



1 - <5(r) 



(C+6) 2 /2 



where the hats indicate that analytic continuation was imposed: 



6+i 



x = - 



1 

6 



and S(r) 



1 — r 
1 + r 



(3f 



Here we may use the momentum conservation in (spatial) directions parallel to the 
Dp-brane and the on-shell conditions, which lead to 



_p\ — t\\ 



(39) 



Notice that the result seems to be asymmetric under the reordering 1 <-> 2 of the 
vertices unless £ = £ 2 as suggested by the analysis of the pole positions in Section I2TT1 
In general £ 7^ £ 2 unless the brane is space-filling (p = 25). We will now focus on the 
symmetric case £ = £ 2 , which leads to a drastic simplification of the result fl37|) . By 
(139]) this then also sets j/c^l = |fcf |. Let us define the kinematic variables 



s 
t 

u 



ki ■ k 2 
2(M') 2 



2(4 



8 - (kty - (kty 



(40) 



which are motivated by the standard conventions (see Appendix[A|) for scattering from 
a stable D-brane [37,44] but should not be regarded as true Mandelstam variables!^*! 
We have crossed s <-» t to obtain a setting which corresponds to string production by 
the brane rather than scattering off the brane. 



10 These variables are not directly linked to 2 — > 2 scattering which is the case for the stable 
D-brane (Appendix For example, even after fixing £ = £2 there are three independent kinematic 
variables in the present case. Hence it cannot be expressed in terms of s and t in contrast to usual 
2-^2 scattering. 
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For equal energies (1371) is apparently singular, because N is continued analytically 
to the value — 2£ where the integral has a branch point. Hence we need to check the 
details of the analytic continuation carefully. We choose a natural branch where the 
continuation is done from positive values of N via the path 



N = 2£e 



lip 



(p = 



71 



(41) 



Then, in particular, r c — > oo along the positive imaginary axis, or more precisely the 
path Re f c = 1/2. The amplitude becomes 



An 




2 sin 2tt£ 



S_ t'2 , t_ rj 

dxx 2 ? ( 1 — x) 2 



dy 



/i-« 2 (l + y)t- 2 -yf(l + y)|-« 2 - 2 



(42) 



where we substituted x 



and y 



After using identities of the gamma 



functions, the result simplifies to 



.4, 



■(2ttA) 



2iw 



2 sinh 2tiuj 
(2ttA) 2 ^ 



\2 2 



1 



B 



s i u 

- + 1, 1 

2 ' 2 



•7TS 

nte 2 



2 sinh 27tcg> 



^ rq- 1 ) , r(§ + i) 
r(^- 2 ) + r(^) 



2 



(43) 



where = c^i = cu 2 = i^. 

Notice that (jl3ll has the expected s- and t-channel poles at s — —2, —4, —6, . . . 
and at t = 2, 0, —2, . . . which arise from the endpoints at r = and r = 1 of the 
integral in (1371) . respectively. It is possible to check that the residues of the first few 
poles are in accord with the general expectations of Section 12.11 The last gamma 
function produces an additional pole at u — 0, while at u = 2 the expression in 
square brackets vanishes which cancels the next pole. The following poles at negative 
u are unphysical since u > by definition. The u = singularity corresponds to 
production of strings which have parallel momenta in the directions perpendicular to 
the Dp-brane, kf = fcf , p < 25. The amplitude is also singular at oj = 0. 



4.1 Kinematic limits 

To conclude our analysis, let us consider the asymptotic behavior of the £ = £ 2 
amplitude (1431) in different kinematic limits. We notice that the kinematic variables 
can be written as 

s = 4-2(£H) 2 -(£ ± ) 2 (l-cos#) 

t = 8-2(^) 2 

u = (k L ) 2 (l -cos6) , (44) 
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xT.l , x 24 



x,'...,x p 

Figure 2: This figure illustrates the kinematic variables: the spatial momenta of the 
emitted strings ki with the angles 6, if). 

where {k^) 2 = (k^) 2 = (k^) 2 and 9 is the angle between k± and k^ (see figure [2]). 
Before going to the explicit asymptotics, we make an important general observation: 
the amplitude in (j4"3|) vanishes exponentially ( ~ e - 27rw or faster) in the high energy 



limit. Recall that one-point amplitude [6] behaves as ~ e" 
since the sums of the arguments of the two beta functions in 
above, 

t + u „ 



This result follows 
are bounded from 



s + u 



2-{k x 
2 - (£H 



, 1 + cos 9 



< 2 



< 2 



(45) 



which excludes the region where the beta functions would grow exponentially for 
u) — > oo. They can have a power- like behavior in u if at least one of the sums is fixed 
for uj — > oo, in which case the behavior ~ e -i™ ar j ses f r0 m the factor 1/ sinh27ro; 
that multiplies the beta functions. The same factor is responsible for the suppression 
of the one-point amplitude. In both cases the exponent contains the total energy of 
the emitted strings. In particular this suggests that in the brane decay the two-string 
emission channel will be as important as the one-string channel. However, we would 
need to study decay into massive string states and the phase space factors before 
drawing conclusions. 

Let us consider the following special limits: 

I (fcH) 2 — > oo with (At - ) 2 fixed, corresponding to string production parallel to the 
decaying brane 

II (Ar - ) 2 — > oo with (P) 2 (and 9) fixed, corresponding to string production per- 
pendicular to the decaying brane 
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Ill (k L ) 2 — > oo with (A;") 2 and ik L ) 2 {l + cos#) (or equivalently k^~ + k^) fixed, 
corresponding to string production perpendicular to the decaying brane such 
that the backreaction to the brane remains finite. 

Notice that all these limits are reasonable, since u 2 = k 2 — A — > oo in each case validat- 
ing the electrostatic approximation. Recall that any subleading terms of the electro- 
static approximation were not included, and hence results for subleading asymptotics 
(which we shall anyhow write down) are not reliable. Also notice that there are 
plenty of poles in the physically meaningful region. We only state the leading overall 
behavior "averaged over the poles" , so the subleading terms may be significant very 
near their poles. 

The limit I means taking s — > — oo with t and u fixed. For u > 2 the factor in 
the square brackets of fj43l) approaches a constant since the latter, s dependent term 
decays, 

A ^ z7re^(27rA)^^r(|-l)r(|-l) 



2sinh27rw T - 2) 

while for < u < 2 the latter term dominates and we find 



(46) 



nre^(27rA) 2 - p /« _ x sin[vr( S + u)/2] /_*\HI ^ (47) 



2sinh27ra> V2 / sin[7rs/2] V 2 
Thus for any u the absolute value of the amplitude vanishes exponentially ~ e~ 27TW , 

2 

as expected. The amplitude has also a very rapidly varying phase ~ _ The 

subleading behavior which arises from the terms in the square brackets is likely to be 
modified by the nonleading terms of the electrostatic limit which we have neglected. 

The limit II corresponds to s, t — > — oo while u — > +oo such that t + u — ► — oo and 
s + u = 4-2(fcH) 2 remains fixed. Then the first term in the square brackets vanishes 
exponentially, and the contribution from the second term leads to 

2sinh2vrcj r (^±m) sin f V 2 J K ' 

Hence the absolute value of the amplitude behaves as ~ e ~ 27TUJ a i so j n this case. 

The limit III is otherwise the same as the limit II except that t + u = 8 — (k ± ) 2 (l + 
cos 9) is fixed. The asymptotics is given by 

_ i7re i ^(2n\) 2w n / s\^-i _ 1 + cosfl (fc") 2 

2 ~ 2sinh27ru; r(^)sinf V 2/ ' 2 > (£l)2 ^ ^ 

i7re^(27rA) 2 ^ vre^ f t\^ -3 1 + cos# < (£H) 2 



2sinh27rcj r -2)sinf V 2 ; 2 " (jfei)2 ' 

That is, in the former case we reproduce (1481) . Again, the absolute value of the 
amplitude vanishes as ~ e~ 2nuj . 
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Let us finally consider string production in an arbitrary direction. We take u — > oo 
with the parametrization 

(P) 2 = fc 2 cos 2 ^ = (4 + cu 2 )cos 2 ^ 

(k^f = k 2 sm 2 ij = (4 + cj 2 )sin 2 ^ , (50) 

where < ip < it/ 2 (such that the above cases I-III correspond to the endpoints, see 
Figure |2j). We find the leading one of the following asymptotics: 



A 2 ~ ^nXy^e-^e-^ 2 ^ 2 ^ 
A 2 ~ ( 2 7rA) 2 ^e" 2 ^e"^ 2 ( cos2,/ ' + 



1 + cos 9 1 — cos i 

l + cos6>\ 2 /l-cos6»^ 2 



-/ , i l-cosfl „ in 2 



1>) 



S i n 2 ^ 



u> sin i/> 



sin 2 1/. 



cos 2 ^ + ^f^sin 2 ^, 



X 



cos 2 ip 



cos if)- 



cos 2 ^ + sin 2 V 



(51) 



The amplitudes vanish very fast ~ exp(— Cuj 2 ) due to the terms with the square 
brackets. It is easy to check that the coefficient C is positive. 
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A Correlators and kinematics 

In this appendix we write the bulk two-point amplitude in terms of the Mandelstam 
variables often used in the literature. Let us start by giving a heuristic argument 
to motivate the definition of the variables. For the moment let us ignore the energy 
transfer from the D-brane to the strings which is possible due to the half S-brane 
deformation. Consider the 2-^2 scattering tachyon (fei) + D-brane (pi) — > tachyon 
(&2) + D-brane {p 2 ), with the associated momenta in brackets. Breaking into parallel 
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and perpendicular components, 



h 


— (4>4~. 


pi 


= (o,-4 


k 2 


= (4')4~. 


P2 


= (0,-4 



(52) 

On-shell conditions and momentum conservation in parallel directions give 



4 


= 4=4 


4 


= 4 


kt\ 


= 1^2 I • 



(53) 

So we obtain the Mandelstam variables 

s = (h + Pl ) 2 = (h + P2 ) 2 = (4) 2 = (4) 2 = \ s 

t = (h - k 2 f = (p! - p 2 ) 2 = (4 X - 4 ± ) 2 = 4 + 4 - 2fci • k 2 = 8 + 2t 

u = (kt - p 2 f = {k 2 - Pl f = (4) 2 + (4" + 4 ± ) 2 = 16 - t - 35 , (54) 

where the variables without tildes match with the conventions of [14] The last line 



of f[54|) is consistent with the relation 

5 + t + u = sum of squared masses = 8 + (4~) 2 + (4~) 2 = 16 — 25 . (55) 
Consider then the bulk two-point tachyon amplitude 

A 2 = J dx°dx lie ^ +k ^l f (T[e ika ' X ' (Wa '* a) \ (56) 



a=l \a=l 

(-z) N r tj du #w a 



jAT /. -N ,, 2 r> / A' 2 

N=l "' " i=l " o=l " \i=l a=l 



+ e / n §? n ^ ( n n ^ dp( -* ) 



The singular self-contractions will be dropped. The first term in the series would give 
the amplitude for scattering from a stable D-brane. It contains the correlator 

Kl = ^JJe^'^-^A . (57) 
Working out the contractions gives 

Kt = exp | - y^{4 • 4'[ ln \ w a ~ w b \ +ln |1 - u> a w fe |] 

a, 6 

+4" • fc^[ln |u> a - w 6 | — In 1 1 — w a w b |]}| . 



11 Notice that has opposite sign in [14]. 
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For the next term in (|56|) we need to evaluate the correlator 

N 2 



e X'°{ti) Y\ P ik a -X'(w a ,w a ) 
\i=l a=l 



We find 



^ = n ~ eitj i 2 n i 1 ~ w ae~ iu \ 2ia ■ k x 



(5f 



(59) 



where £ a = —iuj a . Let us naively generalize the above definitions of the Mandelstam 
variables to the case where the energy transfer is not zero by 



ki ■ k 2 
(*}) a - {ki f 



s 

£-4 



(60) 



where we crossed s «-> t to obtain a setting that corresponds to the production of two 
strings. Then, setting w% = r, wi = 0, we get for the amplitude 



/oo 
N=0 



(61) 



where 



N 



T(N = 0) = T(iV = 1) = / rfrr #i(r) = / drr r s (l - r 2 ) ( ' /2) ~ 2 (62) 
and for N > 2, 

T(N)= ! drrr^l-r 2 )^- 2 ^—^ ! f[ ^ J[ {e* - e*\ 2 J] |1 -re^f* . (63) 

^ ' ^ i=l i<j i 

B Schur polynomial method 

Let us first discuss some properties of the Schur polynomials. They are symmet- 
ric polynomials of N variables x±, . . .xn, denoted by s\(xi,X2, ■ ■ ■ ,%n) where A = 
(Ai, A2, • • • , Xn) is a partition of |A| = Ai + A2 + - ■ - + Atv with Ai > A2 > ■ ■ • > Xn > 0. 
The partition may be identified with a Young diagram where the zth row has Aj boxes. 
The polynomials may be defined by 



det 



SxlXij 



x 



Xi+N-j 



i,j=l,...N 



det 


N-j 








i,j=l,...N 



(64) 



where the denominator is the Vandermonde determinant A(xi, X2, ■ ■ ■ , x^). The def- 
inition may be extended to partitions having more than N nonzero elements Aj (or 
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to Young diagrams having more than N rows) by setting these polynomials to zero. 
We denote by £(A) the index of the last nonzero element (or the height of the corre- 
sponding Young diagram). 

The polynomials have the following well known identities [47] 



A'.r 



nn( i+x ^) = ^2 s *( x i) s y(yi) ( 6s ) 

A 

^sxix^sxiyi) , (66) 



=1 j=l 



nn r 



where A' is the conjugate partition of A, i.e., A^ are the number of boxes in the columns 
of the Young diagram of A. Notice that since s\(xi) vanishes if £(X) > N x , the sum 
over A is restricted to a finite set of partitions with £(X) < N x and Ai < N y in the 
first identity, as expected since the left hand side is a polynomial in X{ and yj. One 
can also write down a generalized binomial theorem 

TV 

Y[(l- Xl Y = Y^c x (Os x ( Xi ) , (67) 

i=l A 

where the coefficients ca(£) are independent of N and polynomials in £. They are 
well known and are related to fl65|) . When £ is a positive integer we can set in fl65l) 
N y = £ and y\ — — 1 = y-i — • ■ ■ — y% whence we see that the coefficient ca(£) equals 
S\'{yi — — 1)- The result may be expressed in various forms: 

ca(o = (-i) |a '^=i)=(-i) iai n K ~ l ~^ +J 

a| a(a-+i,a-_ 1 + 2,... ; a; + 

1 ' (£-l)!---l! 
= C- 1 )'^ II + -i(r^)) , (68) 

?n£A 

where the last form gives the number of terms in the Schur polynomial in terms of 
the "hook length" h m (it also equals the number of the semistandard Young tableaux 
of the form A' and the dimension of the representation of SU(£) characterized by A'). 
Here i(m) (j(m)) is the number of row (column) of the box m in A. Since c\ are 
indeed known to be polynomials in £, the result in the last form readily extends to 
all complex values of £. We may thus write the generalized binomial theorem as 



f[(l - rxif = J2(~ r ) lXl II 7T ( ^ + ' (m) ~ j(m)) Sx ^ ' (69) 

i=l A 

where we rescaled all variables by r. 
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Let us then concentrate on the two-point function Z 2 of (JHJ). By using the formulae 
(I69p we may write the "cross-term" in the integrand as 

N N 

— re-*!* = - re*) 5 (1 - re"*)* = ^(^^^^Mx^s^) , 

i=l i=l A, A 

(70) 

where Xi = e**». As the Schur polynomials are orthogonal with respect to the weight 
function of Z 2 we immediately get 

Z 2 (r,£;N)= ^ |A '|ca(0| 2 = E ^ II ^ + < H ~ J'H) 2 , (71) 

\,£(\)<N X,£(X)<N meA m 

where the constraint £(A) < iV appears since the polynomials sa(xj) vanish for £(X) > 
N. If £ is a positive integer the sum over A is finite because the coefficients ca(£) 
vanish for Ai > £. Moreover, if £ is a negative integer and N > |£| the sum can be 
done explicitly as a special case of the general result found in [24]: 

Z 2 (r,^N) = (l-r 2 )-t 2 ■ f = -1,-2,..., -N . (72) 



C Solutions to the potential problems 

We are interested in solving potential problems in classical electrostatics where the 
conductors are pieces of the arc of the unit circle. The general potential problem has 
as conductors n separate arcs of the unit circle. We find it useful to map the disk 
onto the upper half plane by the conformal map 

w i — ^ z — i/3 = q(w) (73) 

l + w 

with the understanding that we work in the compactified complex plane with the 
oo-point included. Here (3 is a real parameter. 

The analytic properties of the complex electric field allow us to write down the 
solution in a special form [47]. Let us assume for a moment that we have two con- 
ductors, n = 2, which lie at [a, b] and at [c,d\. We can then define the following 
functions 

g(z) = y/z — aV z — &V z — cV z — d 

h(z) = y/z — a\Jb — zy/c — zV d — z , (74) 

where the principal branch of the square root function is used. Then g(z) is ana- 
lytic everywhere except on the conductors where it has branch cuts, whereas h(z) is 
analytic over the conductors, but has branch cuts in the gap regions. 
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Generalization to n conductors is found by adding more terms to g(z) and h(z). 
We choose the branches such that 



h(z) = T ig{z) , (75) 

where the minus sign holds in the upper half plane and the plus sign holds in the 
lower half plane. It is then possible to show that, when exposed to an appropriate 
external (conjugate) electric field E ext (z), a solution for the (conjugate) electric field 
due to the conductors reads 

E c {z)= l -g{z) ! dt ^ E ;*f , (76) 

where the integration is over the conductors. Notice that the solution is bounded 
at the conductor endpoints. It is well defined everywhere except on the conductors, 
where it has discontinuities. The charge density^! is proportional to the discontinuity 
and becomes 

p(z) = —h{z) / dt-— . 77 

7T 2 icond h(t) Z-t 

The 2n endpoints of the conductors are not arbitrary parameters but are subject 
to constraints. In particular, these parameters (the constants a, b, c, and d in the 
h = 2 case) are partially fixed by the asymptotics of the electric field 

E c {z) ~ j (78) 

for z — > oo where N is the total electric charge on the conductors. This gives 

f ^ReEextW , n . 1 

= / dt — , k = 0,...,n — l 

icond hit) 

f ,f ReO) / v 

tiN = dt extv ; . (79) 

icond h{t) 

There are n — 1 parameters left free, which are identified as the potential differences 
between the conducting planes. We will set these to zero and thus require 

Re / dzE(z) = 0, k = 2,...,n, (80) 

where 

E{z) = E c {z) + E cxt {z) (81) 

is the total electric field and the curve connects the first conductor with the kth 
one without intersecting the other conductors, for example. 
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For brevity, we omit the subscript which was used in the main text. 
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For the case of interest to us, the external field is due to n point like charges, 

n 1 

E ext (z)=J2U^— ■ (82) 

z z n 



a=l 



For this configuration the integrals above can be calculated exactly by using contour 
integration methods. We find 



E c (z) 
E(z) 



gjz) 

2 



+ 



glgj r 
9 Z^ Q 



_#(z a )(2:-z a ) #(z*)(z-2:*) 
1 1 



+ 



Hz) 

2tt 



g(z a )(z-z a ) g(z*)(z - z* a ) 
1 1 



^%-z a ){z-zt) 

+ V£ ^ 



53) 



where z a = x a + iy a . The asymptotic conditions imply 







+ 



fc = 0, . . . , n — 1 



1 



_|_ \^a) 



L9{Za) 9{z*)j 



The potential and the total energy of the system are of special interest to us. We 
can integrate the complex electric field in (|83|) to give 



V c (z) 
V(z) 



[Z — Zn 



V(z)+J2^og 

a 

2 Yl ^ a t~ log ( z ~ Za ^ + log 

a 

+ ^{ir* 9(t, ? 5 ° 

iv + 53 e. J log/? 

^ a / 



Z - Z„ 



+ 



p(Za)(*-^a) g(z*)(t - Z*) 



55) 



Notice that this solution has the asymptotics required by our definition ([TBI) (or more 
precisely by its natural generalization to the half plane) 



V(z) 



56) 



so that the constant term of the expansion in 1/z is zero. 

Let us then write down explicitly the solutions for the disk, i.e., apply the transfor- 
mation q(w) of (1731) (with, say (3 = 1). To do so, we should start with the configuration 
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where charges £i> ■ ■ ■ > £n are located at arbitrary points Z\ = q(wi), . . . , z n = q{w n ), 
and an additional charge £ n+i = — N — ^ a £ a is located at z = q(oo). Let us denote 
the corresponding electric field and potential, as given by eqs. ( 1531) and ( 1551) . by E 
and V, respectively. 

The electric field on the disk may be then written as 



E(w) 
EJw) 



E{q{w)) 

E c (q(w)) - 



N 



q(w) - q(oo) 



57) 



and integration gives the potential on the disk 



V(w) = V(q(w)) — lim 

R— >oo 



/ dzE(z)+ \N + )\ogR 

Jm \ a J 



V c (w) = V c (q(™)) + N\og[q(w)-q(oo)] 

fR 



— lim 

R^oo 



dzE{z) -NlogR 



q(oo) 



-N\og-q(l/t) 



t=o 



Here the integrals are independent of w. They were added to achieve the correct 
asymptotics, 



V{w) 



VJw) = -Nlogw + O ( — ] . (89) 

w 



Finally, the potential on the conductors is found as 

U = Re V(w ) , (90) 
where Wo is any point on the conductors, and the "interaction" energy is given by 

£ c -£ = ^£aReKK) • (91) 



The total energy then reads (see (|T9l) above) 

iV 1 



where 



£tot — ~^U + ^Zc-i + £■£ , 



£{ = -^2 Ca6 log(w a - W b ) 

a<b 



(92) 



(93) 



is the energy linked to the self-interactions of the external charges. This trivial term 
was not discussed in the main text since its contribution was embedded in the explicit 
wi,W2 dependence of (jSJ) that was treated separately. 
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